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ABSTRACT 

This paper reports on two XMM-Newton observations of the Seyfert 1 Galaxy H0557- 
385 obtained in 2006, which show the source at an historical low flux state, more 
than a factor of ~10 lower than a previous XMM-Newton look in 2002. The low flux 
spectrum presents a strong Fe Ka line associated to a Compton reflection continuum. 
An additional spectral line around 6.6 keV is required to fit Ka emission from Fe 
XXV. The spectral curvature below 6 keV implies obscuration by neutral gas with a 
column density of ^8xl0 23 cm -2 partially covering the primary emission, which still 
contributes for a few percent of the soft X-ray emission. Absorption by ionised material 
on the line of sight is required to fit the deep trough below 1 keV. The comparison 
of the two spectral states shows that the flux transition is to be ascribed entirely to 
intervening line-of-sight clouds with high column density. 
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1 INTRODUCTION 

A substantial fraction of Active Galactic Nuclei (AGN) 
shows absorption from neutral matter in their X-ray spec- 
tra (Matt et al. 2002 for a review). If the column den- 
\ sity of this gas is higher than 1.5 xlO 24 cm -2 , the level 
of obscuration is such that below 10 keV the X-ray spec- 
trum is dominated by the reflected radiation and the nu- 
cleus is only visible above 10 keV. Transitions between 
Compton-thick (Nh>1-5x10 24 cm -2 ) and Compton-thin 
(Nh<1-5x 10 24 cm -2 ) states have been reported in some 
obscured AGNs observed in the X-rays in the recent years 
(Matt et al. 2003; Guainazzi et al. 2005). The changing look 
phenomenon may be explained in terms of column density 
variations, which are frequently observed in type 2 sources 
(Risaliti et al. 2002) or in terms of the fading of the nuclear 
emission, as proposed for NGC4051 (Guainazzi et al. 1998). 

The Seyfert 1 galaxy H0557-385 (05 h 58 m 02.0 s - 
38°20'05", z=0. 03387) is also known as ESO 055620-3820.2. 
In the X-rays, it was observed by ASCA (Turner et al. 1996), 
BeppoSAX (Quadrelli et al. 2003) and XMM-Newton (Ash- 
ton et al. 2006, hereafter A06). In all these observations, the 
source flux was l-4xl0 -11 ergs cm -2 s _1 and the spectrum 
appeared to be characterised by a strong warm absorber 
below 1 keV. A06 found that a two-phase ionised medium 
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provided a good description of the high resolution data ob- 
tained by the Reflection Grating Spectrometer (RGS). 

XMM-Newton re-observed the source twice in 2006, 
catching it at an extremely low flux state. This paper 
presents the analysis of the spectral transition between the 
XMM-Newton looks in 2002 and 2006. 



2 OBSERVATIONS AND DATA REDUCTION 

H0557-385 was observed 4 times (see Table [TJ by XMM- 
Newton (Jansen et al. 2001). Raw data from all observations 
were reduced with the Science Analysis System (SAS) using 
version 8.0. The EPIC data of 2002 were collected in large 
window mode for the pn camera and in small window mode 
for the MOS1 and MOS2, whereas the 2006 EPIC data were 
all collected in small window. 

Background flares were treated according to the method 
described by Piconcelli et al. 2004 so to maximise the signal- 
to-noise of the data. The spectra of the 2002 observations 
were extracted using pattern 0-4 for the pn and 0-12 for the 
MOS cameras. The 0.3-10 keV count rate of 10 counts/s 
is just below the pile-up threshold in the pn camera. The 
source spectra were extracted from a radius of 40 arcsec, 
yielding a useful exposure of 4 and 6.5 ks in the pn instru- 
ment. No flux variability is found in the light curve and since 
the spectra from the two 2002 exposures are consistent in 
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OBSID 


Date 


pn exp 


^soft 


Fluxo.3_2 


n 

L hard 


Flux 2 _io 




(yyyy-mm-dd) 


(fa) 


(10 — 12 ergs cm -2 s — 1 ) 




(10 — 12 ergs cm -2 s — x ) 



0109130501 2002-04-03 4 1.69±0.08 10.37±0.07 1.60±0.03 36.01±0.05 

0109131001 2002-09-17 6.5 1.53±0.07 10.91±0.06 1.63±0.02 43.30±0.30 

0404260101 2006-08-11 40 2.15±0.13 0.42±0.20 -0.09±0.07 3.06±0.12 

0404260301 2006-11-03 52 2.10±0.15 0.37±0.15 -0.021°;°^ 2.60±0.15 



Table 1. Observation Log of H0557-385 as observed by XMM-Newton. The fluxes have been estimated assuming a soft and a hard X-ray 
power law model (0.3-2 and 2-10 keV respectively). 



flux and spectral shape, they have been co-added in a single 
spectral file. 

The same procedure for the data reduction has been 
applied to the 2006 observations. We note here that due 
to telemetry loss, the observation 0404260101 produced two 
event lists (scheduled and unscheduled), that were subse- 
quently merged in a single event list using the SAS task 
merge. Source and background spectra and light curves 
were then extracted from this merged event list. 

An inspection of the light curves for (both) the 2006 ob- 
servations confirms that the source flux does not vary signif- 
icantly within each of the two exposures, so the analysis will 
be performed on the spectra integrated during each epoch. 

As to the RGS data, for the 2002 observations, we will 
refer to the analysis by A06. For the 2006 exposures, RGS 
data were processed with the task RGSPROC according to 
the standard method proposed in the SAS threads. Due to 
the flux decrease in 2006, the source was not detected in 
the RGS images, therefore no high resolution spectra are 
available for this epoch. 

The source was observed with the Optical Monitor tele- 
scope (OM) with optical and UV filters. Photometry is avail- 
able for the following filters: 5430 (V), 4500 (B), 3440 (U), 
2910 (UVW1), 2310 (UVM2), 2120 (UVW2) A. The data 
have been analysed with the standard pipeline OMICHAIN. 
The source is detected in all the images in each band. From 
the fluxes derived from the observed count-rates, the level 
of variability is lower than 15-20% within the two epochs 
(2002-2006). 



3 SPECTRAL ANALYSIS 

In order to apply the \ 2 statistic the EPIC data were re- 
binned so to oversample the instrumental resolution of a 
factor >3 and to have at least 50 background subtracted 
counts in each spectral bin for the 2002 data and 25 counts 
in the 2006 observations. A Galactic column density of 
4xl0 20 cm" 2 (Dickey & Lockman, 1990) is included in all 
the following spectral fits. Errors are quoted at 90% (confi- 
dence level) for one parameter. Spectral fits are performed 
over the 0.3-10 keV energy range. 

Figure Q] displays the 2002 and 2006 pn data in the full 
EPIC energy band (0.3-10 keV): unmistakably, the source 
undergoes an extraordinary spectral change, with a flux de- 
crease of more than one order of magnitude between 2002 
and 2006 (see Table [TJ. The two epochs are therefore refer- 




Energy (keV) 



Figure 1. EPIC pn spectra of H0557-385 in 2002 and 2006. 



eed to as "high and low flux state" in the following and they 
are analysed separately. 

The final best fit models for the pn data (which are de- 
scribed in the next sections), have been applied to the MOS 
spectra in each spectral state. This confirms that the MOS 
spectra are in excellent agreement with the pn data; for the 
sake of brevity, we quote in the paper only the parameters 
obtained from the pn fits. 

3.1 The high flux state: 2002 data 

The XMM-Newton observations of H0557-385 obtained in 
2002 were published by A06. In the analysis presented herein 
we assume the two warm absorbers model proposed by these 
authors based on RGS data. 

The co-added EPIC pn spectrum of the high flux state 
was fitted in the 2-10 keV band with a power law contin- 
uum and a narrow (zero-width) Gaussian line to account 
for emission in the Fe Ka line at 6.4 keV. The best fit pa- 
rameters are: r=1.65±0.02, E=6.40±0.10 keV, line flux of 
1.60±0.70xl0" 5 photons cm" 2 s" 1 , with a 2-10 keV lumi- 
nosity of ~lxl0 44 ergs/s. The Equivalent Width (EW) of 
the Fe line is 33j^ eV. If the line width is left free to vary the 
measured value is cr=0.20±0.10 keV, with EW=97±33 eV, 
in agreement with the values reported by A06, yielding 
X 2 /d.o.f -284/225. If the Fe line is fitted with a relativis- 
ts profile (DISKLINE, Fabian et al. 1989), the statistic 
worsens 

(x 2 / d -°- f -294/226). The non-zero width of the 
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Table 2. Best fit parameters of the high flux state (pn data) fitted 
with the following model: POWER LAW + NEUTRAL N H +2 
WARM ABSORBERS (ZXIPCF) + 2 GAUSSIANS 



2002 spectrum x 2 /d.o.f.=294/226 



r=1.87±0.01 N H (neutral)=0.07±0.01 
N H (warmi)=1.02±°'3° Log£i=2.16±0.06 
N H (warm 2 )=0.60±0.03 Log£ 2 =0.56±0.05 
Ei=6.40±0.11 Fluxi=1.51±0.71 
E 2 =0.83±0.02 Flux 2 =4.50±1.02 



Line fluxes in units of XlO - 5 ph cm -2 s _1 ; Line energies in 
keV; Column densities in XlO 22 cm -2 ; Ionisation parameter in 
ergs cm s^ 1 . 



Fe line may indicate the presence of additional emission 
components from a higher ionisation state. We tried to fit 
the data with two zero-width Gaussians around 6.4 keV 
and 6.67 keV (mean energy of the Fe XXV triplet). The 
fit statistic is x /d.o.f. =290/224, slightly worse than the 
broad Gaussian fit, but the parameters for the Fe XXV line 
are affected by large uncertainties (E=6.56±0.62 keV and 
Flux=1.34lo ^xlCT 5 photons cm" 2 s" 1 ). The evidence for 
the presence of a complex profile or of two Fe line com- 
ponents in this spectrum is nonetheless very loose, so we 
assume the model with a single narrow Gaussian Fe line to 
compare directly the two spectral states (see Section 4). 

When this model is extrapolated below 2 keV, the data 
present soft X-ray residuals, therefore two warm absorbers 
with ionisation parameters and column densities initially 
fixed to the A06 values have been added in the fit. The model 
zxipcf available in XSPEC as local model was used. Based 
on XSTAR, this model was originally developed to describe 
the effect of a partial covering from ionised matter (Reeves 
et al. 2008). In the present analysis the covering factor is 
fixed to 1 so to mimic a warm absorber with full coverage. 
Beside the multiphase ionised absorber, neutral photoelec- 
tric absorption is required by the data at more than 99.99% 
confidence level, consistently to what is reported in A06, so 
we model it with the XSPEC model ZPHABS obtaining 
Njj=7±1 xlO 20 cm -2 . The final best fitting parameters are 
listed in Table [5] We briefly note here that there is a good 
consistency with the analysis by A06, the small discrepancy 
can be probably attributed to the fact that here only EPIC 
data are considered and to the different models used for fit- 
ting the warm absorber. 

The spectrum still shows positive residuals around 0.8- 
0.9 keV; an additional narrow emission line is included in the 
fit at 0.83±0.03 keV. The improvement in the fit is Ax 2 =33 
for 2 free parameters. The reader is deferred to Section [4] for 
the discussion on the identification of this feature. 



3.2 The low flux state: 2006 data 

As reported in Table [1] the two observations of 2006 were 
performed with a 3 months gap in between, making it nec- 
essary to check for any spectral variation between the two 
spectra. In the following, we will refer to the 2006 spectra as 



Aug06 and Nov06. At first glance, the two spectra are very 
similar in shape and flux intensity. 

The spectral shape in Figure [T] suggests that the drop 
of counts from the high state could be due to absorption 
in the 1-6 keV band. We therefore start by analysing the 
longest spectrum (Nov06) and by constructing a model 
where the 2002 best fit is absorbed by neutral gas. The 
addition of a cold column density to account for the ab- 
sorption (model ZPHABS, already described in the previ- 
ous section) does not provide a good description of the data 
(x 2 /d.o.f-896/233). It seems that the nuclear primary con- 
tinuum is not totally suppressed, but it still contributes to 
the soft X-ray spectrum, therefore a partial covering ab- 
sorber is considered. The partial covering is constructed by 
adding a power law convolved with the ZPHABS and CABS 
components in XSPEC. The latter is necessary to produce 
a correct estimate of the continuum normalisation by tak- 
ing into account Compton scattering by optically thin gas. 
This partial covering model is able to reproduce the spectral 
shape correctly, with column density of 8±0.2xl0 23 cm~ 2 
and covering fraction > 93 percent (see Table [3}. 

The remarkable transition to a low flux state uncov- 
ers a prominent Iron emission line complex (see Figure [2}, 
suggesting that Compton reflection could also contribute to 
the hard X-ray emission. We therefore added a pure Comp- 
ton reflection PEXRAV component (Magdziarz & Zdziarski 
1995) assuming solar abundances and a 30° inclination of the 
reflecting material to the observer line of sight. This com- 
ponent is very significant and it improves the fit statistic by 
A^ 2 =268 for 1 free parameter. The reflection fraction (i.e. 
R=f2/27r) estimated as the ratio of the PEXRAV and the 
direct power law normalisations, is R=0.57±0.04. 

Four zero-width Gaussian lines have been included to 
fit the Fe K complex: Ka and K/3 transitions from neutral 
Fe at 6.4 and 7.06 keV, and Ka lines from Fe XXV-XXVI 
at 6.67 and 6.96 keV. The parameters for the detected Fe 
lines are reported in Table O The Fe XXVI Ka and Fe I K/3 
are not detected in the spectra, so they are not included in 
the table. The upper limit on their fluxes were obtained by 
fixing the peak energy to its laboratory value and they are 
respectively 2.2xl0~ 6 and 2.5xl0" 6 ph cm" 2 s _1 . 

The neutral absorber observed in the high state with 
Nh~7x10 20 cm~ 2 , is not detected in the 2006 data, with 
an upper limit on the column density of ~3xl0 19 cm -2 . 

The effect of the soft X-ray absorbers is still very evident 
in the data, so the same warm absorbers used for the high 
state are included in the model. An additional Gaussian line 
is required to fit residuals around 0.9 keV and it has been 
included in the fit (Ax 2 =83 for 2 free parameters). 

This model provides a satisfactory fit for Nov06 and it 
was then applied to Aug06. The best fitting parameters for 
both spectra are reported in Table [3] The two spectra are in 
very good agreement, therefore we can safely assume that 
the source was not affected by significant spectral change 
between August and November 2006. 

The drastic change in the X-ray emission could also 
be explained if the primary emission simply drops, letting 
only the reflected spectrum be observable (Matt et al. 2003). 
This hypothesis is tested by fitting the November data with 
a PEXRAV component to account for the hard X-ray spec- 
trum, and a single power law to model the soft X-ray excess. 
The warm absorber components and all the Gaussian lines 
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Table 3. Best fit parameters of the low flux state (pn data) fit- 
ted with the following model: POWER LAW +2 WARM AB- 
SORBERS (ZXIPCF) + PEXRAV + PARTIAL COVERING 
(CABS*ZPHABS*PLAW) + 3 GAUSSIANS 



AUGUST 2006 



X 2 /d.o.f.=249/217 



r _i 94+0.10 
N H (pc)=81±2 
N H (warmi)=1.02± o -jjj 
N H (warm 2 )=0.10±°;* 



E 2 =6.66±0.06 
E 3 =0.89±0.01 



R=0.73±0.01 
C/>0.93 
Log6=2.28l°-^ 
Loge 2 =0.28±g;g 
Fluxi=1.08±0.17 
Flux 2 =0.26±0.15 
Flux 3 =1.52±0.36 



NOVEMBER 2006 



X 2 /d.o.f.=263/232 



r _l nc+0.02 

i — i.ao_ 19 

N H (pc)=80±2 

N H (war mi )=1.25l°;J| 

N £r (warm 2 )=0.10±0.01 
El _ 6 41 +0.02 

E 2 =6.62±0.03 
E 3 =0.90±0.01 



R=0.57±0.04 



Log£i=2.15±0.07 
Log€ 2 =0.28l°;°| 
Fluxi=1.40±0.15 
Flux 2 =0.42±0.13 
Flux 3 =1.60±0.28 



Line fluxes in units of X 10 ph cm 



Line energies in 



keV; Column densities in xlO cm ; Ionisation parameter in 
ergs cm s~ . 




Energy (keV) 

Figure 2. EPIC pn spectra of the low flux state in the Fe K 
band. The dotted lines mark the position of the Ka transitions 
for Fe I and Fe XXV (6.4 and 6.67 keV). 



reported in Tableware included in the fit. This model yields 
an acceptable fit statistic (x /d.o.f.=274/233), but the in- 
trinsic photon index of the reflected component is unrealis- 
tic (r=1.13±0.03). If the PEXRAV photon index is fixed to 
the high state intrinsic power law (F=1.9), the \ 2 worsens 
considerably, being x 2 /d.o.f.=458/234. These results allow 
us to exclude that H0557-385 was reflection dominated in 
the 2006 low state, therefore this interpretation will not be 
discussed in the remainder of the paper. 



4 DISCUSSION 

The spectral analysis in the previous section reveals a pe- 
culiar behaviour for this Seyfert 1 galaxy. From a general 
overview of Tables [2] and [3] it is clear that most of the fit- 
ting parameters do not vary between the 2002 and the 2006 
spectra. The photon index of the power law and the warm 
absorbers column density and ionisation state are remark- 
ably consistent within the two epochs. Notably, the 2002 
data require a simpler model (Table [2} compared to the one 
found for the 2006 data in Table [3] Nonetheless, for a more 
direct comparison of the two states, the low state model is 
applied to the high state spectrum, as it follows. 

The normalisation of the power law at 1 keV 
is around 1.45±0.05x 10 -2 ph cm -2 s _1 and 
1.02±0.05xl0" 2 ph cm" 2 s" 1 for the high and low 
state, respectively. These numbers imply that the intrinsic 
X-ray photons flux remains visible in both states with a 
variation lower than 50 percent between 2002 and 2006. As 
described in Section 13.21 a fraction higher than 93 percent 
of the intrinsic flux is absorbed by circumnuclear material 
with high column density (~ 8xl0 23 cm -2 ) during the low 
state. The high state data do not show evidence for partially 
covering gas. When applying the low state model to these 
data and after fixing the column density to 8xl0 23 enrt 2 
(i.e. the low state value), the upper limit on the covering 
fraction is about 5 percent. The emerging picture is that 
of an active nucleus covered by intervening optically thick 
cold gas in 2006, which may have been present during the 
high state epoch with a lower column and/or with sparse 
coverage, as tested by applying the partial covering model 
to the 2002 data. 

The flux of the narrow Fe Ka line at high and low state 
are in good agreement (Tables [2] and [3}. The absence of 
a Fe Ka broad relativistic component and the lack of Fe K 
variability seem to point to a scenario in which the Compton 
reflection occurs on the parsec-scale torus-shaped material 
expected to be present in Seyfert galaxies (Antonucci 1993). 
The Fe line intensity measured with respect to the reflec- 
tion component is EW=460t2ig eV in the high state and 
EW=510±55 eV in the low state. Both values are consider- 
ably lower than the theoretical 1 keV Fe line expected from 
reflection in a Compton-thick torus (Ghisellini et al. 1994). 
Alternatively, a substantial fraction of the Fe line could be 
emitted by transmission in the obscuring clouds. Anyhow, it 
is very difficult to separate Fe line emission from the torus 
and from the clouds without a clear picture of the AGN 
geometry. 

We can speculate on the distance of the clouds by esti- 
mating their velocity from the FWHM of the Fe K line at low 
state. The upper limit on the line width is 45 eV, resulting in 
a velocity of 5000 km s -1 . The source is known to have sev- 
eral components of the Ha line (A06 and reference therein) 
with FWHM of 1035, 2772 and 11000 km s~\ suggesting a 
location of the obscuring clouds consistent with that of the 
two lower velocity components of the Broad Line Region. 
Assuming that in 2006 the source remained at low state for 
3 months, this gives a minimum duration of the occultation 
event. The clouds can then be located at a distance higher 
than vxt=5000 km s _1 x3 months i.e. ~3.9xl0 15 cm. 

Determining more detailed properties of the system 
such as the precise time scale for the transit of the cloud 
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or a characteristic size, appears very difficult since H0557- 
385 was observed very sparsely by X-ray observatories. The 
closest look to the XMM-Newton ones, is an ASCA obser- 
vation of 1995 (Turner et al. 1996) which shows the source 
at high state. In the attempt to find additional X-ray data, 
H0557-385 has been searched for in the XMM-Newton Slew 
Survey images (Saxton et al. 2008); luckily, one recent slew 
(25/02/2008) passed through the coordinates corresponding 
to the AGN (Saxton, private communication). H0557-385 
was not detected, but a reliable upper limit on the flux 
could be derived: <1.5xl0~ 12 ergs cm -2 s -1 for the 0.2- 
2 keV band and <lxl0~ n ergs cm" 2 s" 1 in the 2-12 keV 
band. Arguably, this demonstrates that after ~2 years, the 
source was still at low state. Whether it has remained in 
this state since 2006 through February 2008 is, obviously, 
impossible to tell. 

From the analysis of the UV data (Section 2) it is found 
that the source of UV photons remains essentially unaffected 
by the X-ray change. This implies that the size of the X-ray 
source must be very compact compared to that of the UV 
source. 

Highly obscured AGN constitute an excellent labo- 
ratory to study the physical properties of the circumnu- 
clear gas via soft X-ray spectroscopy of emission lines pro- 
duced in the interaction between such gas and the pri- 
mary nuclear emission (Guainazzi & Bianchi 2007). Un- 
fortunately, the high resolution cameras did not provide 
a detection for H0557-385 at low state and so the source 
does not have the wealth of additional spectral informa- 
tion that was available for other Seyfert 1 caught at low 
state, e.g. Mrk 335 (Longinotti et al. 2008), NGC 4051 
(Pounds et al. 2004). Nonetheless, the 2006 EPIC data 
do show an emission line around 0.9 keV, which is con- 
sistent with emission from the NelX triplet. We have 
searched for this component in the RGS spectrum of 
2002. Assuming the best fit model in Table [2] the line 
is detected with Flux=3.4lt2 34 x 10~ 5 ph cm -2 s" 1 and 
E=0.916±0.017 keV. 

In addition to the NelX line, the 6.6 keV emission line 
identified as Fe XXV in the low state spectrum could also 
be a signature of photoionisation. However, the statistical 
quality of these measurements prevent us from drawing any 
firm conclusion on the origin of these spectral features. 

Detecting soft X-ray emission (both lines and contin- 
uum) may lead to interpret H0557-385 at low state within 
the typical Seyfert 2 scenario (Bianchi et al. 2006; Guainazzi 
& Bianchi 2007). Remarkably, the presence of the same 
two phase-warm absorber in both states plays a key role 
in ruling out this scenario. It provides a very robust argu- 
ment against the "Seyfert 2/scattering" scenario and con- 
sequently, in favour of the partial covering. The fact that 
the warm gas is still visible through the absorption at low 
state, implies that the soft X-ray emission does not come 
from an extended region as the NLR in obscured AGN (it 
would be very unlikely to have a warm absorber at such large 
scale). Accordingly, the soft X-ray emission must have been 
only partially reduced by the obscuring event, as also sug- 
gested by the comparison of the power law normalisations, 
and therefore it is still observable at low state. 

The fact that the Seyfert 1 nucleus remains visible 
at low state makes the spectral transition in H0557-385 a 
very interesting event, which resembles the occultation phe- 



nomenon seen in the Seyfert 2 galaxy NGC1365 (Risaliti et 
al. 2005, 2007), but it is the presence of the warm absorber 
what makes H0557-385 a unique case where the partial cov- 
ering scenario is indeed very robust and can be challenged 
only if other spectral states were observed. 
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